Macroscopic quantum self-trapping in a Bose-Josephson junction with fermions. 
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We study the macroscopic quantum self-trapping effect in a mixture of Bose-Einstein condensate 
and a large number of quantum degenerate fermions, trapped in a double-well potential with a 
variable separation between the wells. The large number of fermions localized in each well form 
quasi-static impurities that affect the dynamics of the bosonic cloud. Our semi-analytical analysis 
based on a mean-field model shows that main features of macroscopic quantum self trapping in 
a pure bosonic system are radically modified by the influence of fermions, with both the onset 
of self-trapping and properties of the self-trapped state depending on the fermion concentration 
as well as on the type of inter-species interaction. Remarkably, repulsive inter-species interaction 
leads to population inversion of the bosonic energy levels in the trapping potential and hence to 
the inversion of symmetry of the macroscopic wavefunction. Realistic physical estimates are given 
based on experimental parameters for a 40 K- 87 Rb system. 

PACS numbers: 03.75.Lm, 67.85. Pq, 03.75.Hh 



I. INTRODUCTION 

Self-trapping phenomena are among the most dramatic 
effects of atomic interactions [lj in the systems of quan- 
tum degenerate gases. The so-called macroscopic quan- 
tum self-trapping (MQST) effect @ manifests itself as 
localization of most of the particles in the system in a 
particular region in space. The MQST and related ef- 
fects in effectively purely bosonic systems have been ex- 
tensively analyzed in different physical contexts, from the 
Josephson effect in superconductors [3( and the study of 
surperfluid He 4 [4] to the alkali Bose-Einstein conden- 
sates @. 

The MQST effect in a so-called Bose Josephson junc- 
tion, i.e. a Bose-Einstein condensate loaded into a 
double-wellpotential, has been extensively studied the- 
oretically [g, 0] and observed experimentally [§]. Its 
appearance is linked to a nonlinearity-induced phase- 
locking between the nonlinear eigenstates of the system, 
which leads to the formation of new stationary states in 
the trapping potential that have no linear counterpart 
0, H|. These states are characterized by a population 
imbalance between the wells of the trapping potential, 
which becomes more pronounced with growing nonlinear- 
ity (see Fig. Q}. Moreover, it turns out that the MQST 
effect plays an important role in the dynamics of conden- 
sates in periodic potentials, leading to the formation of 
self-trapped or truncated gap states [13, HH • 

The experimental achievements in the production of 
mixtures of BECs and degenerate Fermi gases with inter- 
species interaction strength controlled by Feshbach reso- 
nances [TU, 0, [3, EH , drew attention to the studies of 
the Bose-Fermi systems in the context of novel quantum 
phases and lattice dynamics. Properties of the nonlin- 
ear spatially localized states, Bose-Fermi solitons, have 
also been studied both in harmonic 

nana 

Ull and pe- 
riodic trapping potentials fl9| . In contrast, this paper 



presents a theory of the self-trapping effect in the basic 
model of a Bose-Fermi mixture loaded into a double-well 
potential. As demonstrated by the previous research into 
purely bosonic systems [13 HH , such an effect will also 
have consequences for the behaviour of the quantum de- 
generate mixture in periodic potentials. 




FIG. f: Schematics of the atomic density distribution for 
fermionic impurities (dashed line) and bosonic condensate 
(solid) in a macroscopic self-trapped stationary state formed 
in a one-dimensional double- well trapping potential (dotted) . 

Our aim here is to describe the effect of degenerate 
fermions on the self-trapping behavior of ultracold bosons 
in a quasi-one dimensional symmetric double- well pqten- 
tial [8[ , in the limit of a large number of fermions [201 ] . In 
this limit, the ultracold fermions are not likely to engage 
in a collective motion and, in a presence of a confining 
potential, form a quasi-static inhomogeneous distribution 
of localized impurities in the bosonic cloud (see Fig. []}. 
In the case of a double-well potential we develop ana- 
lytical treatment of the system, based on a quasi-static 
density approximation for the fermionic component and 
a coupled- mode theory for the bosonic component. We 



2 



consider both attractive and repulsive interactions be- 
tween bosons and fcrmions, which can be modified by 
means of Feshbach resonances d [H, El (3 [H, H , in 
particular considering the system of ultracold 40 K- 87 Rb 
atoms. We analyze the effect of a quasi-static distribu- 
tion of fermions on the formation of MQST states in the 
interacting bosonic component of the system and show 
the difference of the physics between the attraction or 
repulsion between fermions and bosons. The two types 
of interaction between the atomic species can be achieved 
experimentally [l2| , which offers a unique opportunity to 
explore different MQST phases in the system. 

The paper is organised as follows: Section II introduces 
the second quantized and the effective one-dimensional 
mean-field model for the quantum degenerate Bose- Fermi 
mixture trapped in a one-dimensional double-well poten- 
tial. Section III presents our basic approximation of the 
quasi-static density distribution for fermions and deriva- 
tion of the coupled mode theory for the macroscopic BEC 
wavefunction. Section IV presents the analysis of the self- 
trapping regimes based on the coupled-mode theory and 
variational ansatz for the BEC wavefunction in the two 
main cases of the inter-species attraction and repulsion. 
Section V concludes and summarizes our work. The de- 
tails of the variational calculations are given in the Ap- 
pendix. 



II. THE MEAN-FIELD MODEL 

We consider a Bose-Fermi mixture trapped in a one di- 
mensional symmetric double-well potential. The Hamil- 
tonian of the system can be written in second quantized 
form in terms of single particle operators for each of the 
particle species and the inter-species interaction between 
bosons and fermions: 

n= / d J + 9 -t^Maa 

where the bosons and fermions in the system are repre- 
sented by their corresponding field operators: i/^ j (cre- 
ation) and ip b j (annihilation); tp^ = p)^{r,t). The field 
operators obey the usual commutation (bosons) and anti- 
commutation (fermions) algebra. 

The longitudinal (^-component) single-particle part of 
the Hamiltonian (TT|) includes a one-dimensional double- 
well trapping potential: H^ x = p 2 /(2m^) + m^ui 2 (\x\ — 
xo) 2 /2. where the trapping frequencies and the masses 
for each of the species are denoted by uj b j and m b j, 
respectively, and the distance between the minima of the 
double- well is given by 2xq. 

The interaction strengths g bb and gbf in Eq. |T]) rep- 
resent the scattering between bosons and boson-fermion, 



respectively. These coefficients are considered to be con- 
stant in the relevant region near their resonances and 
approximately proportional to their s-wavc scattering 
length. The on-site fermion-fcrmion interaction is Pauli 
suppressed. In the following sections we discuss both he 
inter-species repulsion and attraction, and show that the 
sign of the interaction leads to very different physics. 

Following the treatment given in (l9| . we derive the 
equations of motions from the Hamiltonian ([1]) by using 
the Green's function method, and perform a mean-field 
approximation, thereby defining the condensate wave- 
function ip b w (V'fc)' an( i ea ch of the wavefunctions cor- 
responding to the fermions, pfj ~ yPf 1 )- We consider 
the quasi-one dimensional limit of a strongly elongated 
trap with the transverse trapping frequency luj_, which 
leads to rescaling of the one-dimensional interaction co- 
efficients gbb — 2hio±a bb and gbf — 2huj±a b f. Thus, we 
arrive to a set of Nf + 1 coupled mean-field equations, 
Nf equations corresponding to each of the fermions and 
one for the BEC wavefunction [l9j : 

ihdti>] = (H f + ub f p°b)^, (i) 

m t pj° b = (Hb + UbbP°b+UbfP f )Pjl (2) 

where, p° = a J2n=i IV 7 / | 2 j a is a srze parameter of 
the cloud, and p b = \ip b \ 2 - The number of bosons 
(fermions) in the system is given by: Nb.f — J R dxp b ^ . 
By introducing the scaling units of time, 2/uob, length, 
yjh/ '(rribLUb), and energy, hujb, the model can be re- 
written in the following dimensionless form: 

\^-d t r f + d 2 x r f - v f r f + «i# 6 iv? = o, (3) 

idtiPl + dliPl - V b ^ b + a\ip° b \ 2 ipl + u lPf iPl = 0, (4) 

where the bosonic wavefunction and fermionic density 
are rescaled as ip b — > ip b [fvuj b / '(2|ubb|)] 1 / 2 and — > 
p° f fkvf/(2\ubb\), respectively, u = -u b fmf/(\u bb \m b ), 
ui = -u b fujf /(\u bb \u b ), a = -sgn(u bb ). The double- 
well potential for bosons or fermions, = —m^K^(\x\ — 
Xo) 2 /m b , is parametrized by the ratio of the trapping 
strengths, = (m^u 2 ) / (rn b u> 2 ) that determines the spa- 
tial scale of the effective trapping potentials experienced 
by the different species of atoms. 

III. COUPLED-MODE THEORY 

Numerical solutions of the model equations J3]) are 
trackable when the number of fermions is small 11911 . 
However, the current state-of-the art experiments [201 ] 
suggest that the system should be considered in the limit 
of large number of atoms, both for bosonic and fermionic 
component, and the number of degenerate fermions can 
be of the same order of magnitude as the number of 
bosons. We consider the density of fermions, p", tak- 
ing its stationary limit near the condensation tempera- 
ture of the bosons, with the dynamics of the fermionic 
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cloud occuring on a much larger time scale than that 
of the bosonic cloud. In this limit, the fermions can be 
described in a first approximation by a filled Fermi sea, 
where the Fermi points ±kp are invariant upon the inter- 
action, and Luttinger's theorem holds for the fermionic 
component [13, . We assume that this quasistatic dis- 
tribution of the fermionic density takes the form: 



P°f 



p° f (x,0) 



Xq 



(•5) 



where r/& = Nf/Nb is the fraction of fermions in the 
system, which corresponds to the fermions spatially lo- 
calized at the minimum of each well (as shown in Fig. 
[T]). The size of the localization region is governed by the 
parameter a ss \Jmbl '( m / K /)j with larger values of a 
corresponding to greater spatial extent of the fermionic 
cloud. The choice of density function also ensures 
that the number of fermions in each well is conserved. 

The assumption of a quasistatic fermionic density dis- 
tribution introduced above neglects the dynamics of indi- 
vidual atoms given by (O , as well as the feedback of the 
bosons on the fermions, while imposing a boundary con- 
dition on the total density. However, this approximation 
allows for analytical treatment of the model equation (0]) 
governing the behaviour of bosons. 

In accordance with the coupled-mode theory for Bose- 
Einstein condensates trapped in a double-well potential, 
developed in 0], we assume that the main contribution 
to the bosonic condensate wavefunction comes from the 
two lowest nonlinear modes of the double- well potential, 
corresponding to the symmetric ground state, $>o(x), and 
anti-symmetric first excited state, <&i(x): 



where: C i = L ^l^jdx / (nonx) , and tt = 2(m - u ) + 
1u(C\n\ — CqUq). This system of non-liner coupled equa- 
tions describes the dynamical population exchange be- 
tween the lowest energy states of the condensate cloud 
in a double- well potential, and provides an alternative to 
a coupled mode description in the basis of ground states 
of the two separate wells. The advantage of the model 
is in its ability to treat the system for any separation 
between the wells, providing an accurate picture of the 
population dynamics in the situation when the tunneling 
between the wells is strong. As shown in Q, the dynami- 
cal system §8§ admits phase-locked solutions correspond- 
ing to macroscopically self trapped states, characterised 
by the arrest of tunneling and formation of a stationary 
atomic density distribution which is nonzero in each well 
but strongly unbalanced. 

In order to analyze possible regimes of the self-trapping 
in the presence of localized fermionic impurities, we will 
set the physical parameters of the system to those of a 
40 K- 87 Rb quantum degenerate mixture @, GJ, 0, GJ, EH • 
We study the system where the rescaled one-dimensional 
s-wave inter-species interaction parameter is attractive, 
Ubf ~ —234 ao for 40 K - 87 Rb and the analogous repulsive 
case, with Ubf ~ +234 ao, where ao is the Bohr radius. 
The boson-boson scattering parameter is Ubb ~ 98.98 ao 
[2(il ] . The double- well angular frequency for 87 Rb is u>b ~ 
2tt x 263 Hz (see [8j) and u>f = Uby/mb]JjnJ^7) ~ ^b/a. 

In the 



and u>f = ujb yJmb/(mfKf] 
In a realistic experimental situation « 1 \24 _ 
following section, we examine the difference in the self- 
trapping scenario for bosons caused by the different types 
of inter-species interactions. 



IV. SELF-TRAPPING REGIMES 



$ = J2 *j( x ) B M eM-W + o-C jnj t)/^nJ, (6) 



3=0 



where Bj (t) is the time-dependent amplitude of the rele- 
vant state (j = 0, 1), [ij is its energy, rij = J R dx and 
Cj = J* R dx $j/nj. The nonlinear modes $ j obey the 
following stationary equations derived from Eq. (0J: 



dx 2 



+ H j $j-(\x\-x ) 2 $ j + a<f> 3 j +u lP f (x)<f> j =0, (7) 



The equation ([7]) is equivalent to a standard mean field 
Gross-Pitaevskii equation for a BEC with a position de- 
pendent potential term modified by the nonlinear inter- 
action with the fermionic cloud. 

Substituting the ansatz ^ into Eq. ^ and using Eq. 
([7]) , we recover the system of two coupled equations for 
the mode amplitudes 0]: 

[ liT = ° C o\Bo\ 2 B + o-Coi (2|B 1 | 2 B + So^i 2 e~ im ) , 
i^§i = aC 1 \B 1 \ 2 B 1 +aC 01 (2|S | 2 Si + Br^e" int ) ,(8) 



In order to determine the parameter space 
where MQST occurs, we use the ansatz, Bj(t) = 
\J rij (t) exp [— i<f>j (t)], which allows us to rewrite the 
system © in terms of the the difference in relative 
populations of the modes A = n± — no and the relative 
phase shift, ip — 2(</>o — 4>i) ~ £lt [7j: 

dA 2 o . 

— =crG i(n - A )siinp, 
dt 

^- = -6 + a(C + d)A - 2aC* i(2 + cos p)A, (9) 
dt 

where S = 2(/ii — fig) + a[(n — 2no)Co — (n~ 2rii)Ci] , and 
n = no + ni = coiist. The self-trapped states correspond 
to the regime where the relative phase shift is fixed to 
an integer of 27r, ip — 2irm. The difference in the popu- 
lation for the MQST states is given by the self-trapping 
parameter: 



^Co + d-ecoi)' 



-n < A01 < n. 



(10) 



The limiting cases correspond to the situation when only 
symmetric or antisymmetric mode is populated, resulting 
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in the equal number of atoms in each well. The ultimate 
self-trapped state, when all the atoms are localized in one 
well, corresponds to A<ji = 0. 

The value of the self-trapping parameter A i is 
strongly influenced by the mode coupling strengths CV,- 
and the effective energies fij of the nonlinear modes. In 
our model these can be determined semi-analytically by 
means of a variational approach. This approach employs 
an ansatz for the macroscopic wavefunction of the con- 
densate, $, in the form of a linear combination of the 
symmetric $o, and anti-symmetric $1 eigenstates of the 
double-well potential: &(x, Xq) = $o(x,xq) + &i(x, Xq), 
where: 

(_(x-x Q ) 2 _(* + x ) 2 \ 

e "^S7~ ± e "^Tj , (11) 

and the amplitudes, Aj, and widths, dj, are the varia- 
tional parameters to be determined for each fij . The form 
of this ansatz reflects the existence of the first two lowest 
energy eigenstates, with symmetric and anti-symmetric 
spatial profiles and maxima localized in the minima of 
the double-well potential. The details of the calculations 
are presented in the Appendix. Hereafter we present the 
results of the variational analysis. 




FIG. 2: The difference in energies, pi — (io, as a function of 
the double- well separation xo- Inset: variation of the widths 
of the condensate wavefunctions, oo,i, vs. xo- Parameters 
are: r/f, = (solid line), 0.2 (dotted), 0.625 (dashed), and 
a = 0.77. 



A. Inter-species attraction 

We begin by analyzing the difference in energies, \x\ — 
Ho, as a function of the double- well separation Xq (see 
Fig®, in the case of «&/ < 0, i.e. attraction between the 
fermionic impurities and BEC. In this case the growth in 
the concentration of fermions, rn, leads to the growing 
energy splitting between anti-symmetric and symmetric 
state, as compared to the pure BEC. The growth of the 



energy difference leads to the stronger localization of the 
condensate wavefunction in each well (see the inset in 
FigJ^J), which in turn results in the suppression of tun- 
neling and the onset of the MQST at smaller well sepa- 
rations, compared to a pure BEC case. 

The behaviour of the population imbalance, Aoi(xo), 
is shown in Fig[3J In the region marked there is no self- 
trapping effect. In contrast, for region I, self-trapping 
occurs at smaller separation than the pure bosonic sys- 
tem, the latter restricted to the region 77. As we in- 
crease the separation the system moves deeper into the 
self-trapped regime, where the macroscopic wavefunction 
is completely localized in one of the two wells, and both 
the symmetric and antisymmetric mode are equally pop- 
ulated. Spatial profiles of the BEC wavefunction for the 
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FIG. 3: The macroscopic quantum self-trapping parameter 
Aoi/n as a function of the well separation, xo, for attractive 
inter-species interaction. Parameters are: r/b = (solid line), 
0.2 (dotted), 0.625 (dashed), and a = 0.77. The roman nu- 
merals stand for regions with different symmetry properties 
of the BEC wavefunction (see text). 

different regions in the parameter space are shown in 
FiglU We can see that, although the attractive inter- 
species interaction leads to fermion induced narrowing of 
the wavefunction peaks, the effect is weak even for large 
concentrations of fermions. Nevertheless, the difference 
in the strength of the coupling between the nonlinear 
modes of the system produced by the small variations in 
the spatial localization may lead to a significant fermion- 
induced extension of the self-trapping regime. 

The data from the experiments on the Bose-Josephson 
junction [f| provide a reference for experimentally feasi- 
ble well separation. Provided that the Bose-Fermi mix- 
ture can be achieved with a strong inter-species interac- 
tion (i.e. either a wide spread of fermionic impurities, 
given by a, or \ubf\/ao ^S> 1), the detection of the MSQT 
regime in the mixture is possible within the range of xq 
to 6.65 and 10.15, which corresponds to the experimen- 
tally achieved range of well separations from 4.4 /iin to 
6.7 /im. We note that the strength of the nonlinear in- 
teraction between the BEC and fermionic impurities can 
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be adjusted by modifying either the inter-species scatter- 
ing length or the trapping frequencies for the two species 
(u)b ~ aw/). 
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FIG. 4: The macroscopic wavefunction, $(x,xo), for tfb = 
(solid line) and 0.625 (dashed) with attractive inter-species in- 
teraction and xo corresponding to the marked points in Fig[3] 



exist in all regions except for the region 0. The onset 
of the MQST effect in a Bose-Fermi mixture occurs at 
larger well separations compared to that in a pure BEC 
system (see region J). This is due to the fact that the 
presence of repulsive fermionic impurities leads to the 
effective broadening of the BEC wavefunction [see Fig. 
IH1 (inset) and Fig. [7] (a)], which extends the regime of 
strong inter- well tunneling to larger well separations. At 
the well separation corresponding to the point (b) all 
atoms are localized in one of the well, like in the pure 
BEC case. Beyond this point, in region II, new behav- 
ior emerges, where owing to the population inversion be- 
tween the ground and the first excited state, the MQST 
states with nonzero population in both wells emerge once 
again. The bosonic wavefunction corresponding to the 



B. Inter-species repulsion 

The energy splitting between the ground and the first 
excited state of BEC in the double-well potential in the 
case of repulsive inter-species interaction, Uft, > 0, has a 
very different behavior compared to the attractive case. 
In the limit of small separations, the symmetric state 
is the lowest energy state in both pure BEC and Bosc- 
Fermi mixture. In the pure BEC case the energy split- 
ting monotonically tends to zero resulting in the energy 
degeneracy at large separations. However, the depen- 
dence of the energy splitting on the well separation in a 
Bose-Fermi mixture is non-monotonic, with an extremum 
value different from zero, see FigO For a certain range 
of well separations, namely in the region where the differ- 
ence /ii — no changes sign, the ground state of the bosonic 
component of the mixture becomes antisymmetric. 




-6\l . . . , . , . , . I 

1 2 3 4 5 

X 

FIG. 5: The difference in energies, pi — (io, as a function of 
the double- well separation xo for repulsive inter-species inter- 
action. Inset: the variation of the widths, ao,i vs. xo. Param- 
eters are: r/{, = (solid line), 0.2 (dotted), 0.625 (dashed), 
and a = 12.2. 

The dependence of the self-trapping parameter Acu on 
the well separation is shown in Fig. [5] The MQST states 
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FIG. 6: The macroscopic quantum self-trapping parameter 
Aoi/n as a function of the well separation, xo, for repulsive 
inter-species interaction. Parameters are: r/i, = (solid line), 
0.2 (dotted), 0.625 (dashed), and a = 12.2. The roman nu- 
merals stand for regions with different symmetry properties 
of the BEC wavefunctions (see text). 

point (c) in region 77 of Fig. [6] are shown in figure [7] 
It can be seen that the symmetry of the MQST state in 
this region is inversed compared to the pure bosonic case, 
namely the peaks of the wavefunction in the two wells 
are in-phase. This is in a sharp contrast to the "con- 
ventional" self-trapped out-of-phase states, which occur 
both in purely bosonic BEC and in the Bose-Fermi mix- 
ture [cf. Figs. [Ha) and|Hc)]. The fermion induced sym- 
metry breaking of the macroscopic wavefunction is due 
to the fermion-induced population inversion between the 
symmetric and the anti-symmetric mode, n± < no- Pop- 
ulation imbalance vanishes asymptotically at large well 
separations in region II, whereby the bosons once again 
occupy only one of the two wells. 

The highly nontrivial nature of the self-trapping in the 
case of repulsive interaction between the BEC and lo- 
calized fermionic impurities can potentially be explored 
in an experiment with a double- well potential, similar to 
that presented in [f|. In an experiment, however, one 
would measure the difference in atom numbers between 
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FIG. 7: The macroscopic wavefunction, $(2;, £0), for rj b = 
(solid line) and 0.625 (dashed) with repulsive inter-species in- 
teraction and xo corresponding to the marked points in Fig[6] 



the BECs occupying two different wells of the potential, 
AiV, rather than the populations of two nonlinear modes 
$o,i- Possible results of such a measurement are shown 
in Fig. [5Ja). It can be seen that, in contrast to the case 
of a pure BEC (solid line), the population imbalance, 
AN/NT t, where Nj~ot is the total number of BEC atoms, 
is a nonmonotonic function of well separation [dashed line 
in Fig. [5fa)] approximately given by: 



AN 

N T ot 



1 + A§! 



erf(y) 



(1 -cxp(-2y)) 1 /2' 



v = ^ (12) 



At large concentration of fermions the population im- 
balance exhibits a sharp minimum at the certain value 
well of separation [see inset in Fig. [SJb)]. The minimum 
value of the imbalance tends to zero as rjt — > 1 which 
means that, at large well separations, the self-trapping 
effect not only becomes less pronounced with the growth 
of the fermionic concentration, but may disappear en- 
tirely. We should note here that our assumption of the 
static fermionic density localized around the two minima 
of the potential wells becomes less reliable for r/b > 1. 



V. CONCLUSIONS 

We have analyzed the formation of self-trapped states 
in a BEC cloud mixed with degenerate fermions and con- 
fined in a one-dimensional double- well potential. Our an- 
alytical approach is reliable in the limit of large number of 
fermions and is based on the quasi-static approximation 
for the fermionic density, where the fermions are consid- 
ered to form large impurities spatially localized around 
the minima of the two potential wells. The properties 
of the bosonic macroscopic wavefunction are analyzed by 
means of variational method, which allows us to compre- 
hensively describe the spatial properties and symmetry 
of the self-trapped state. The self-trapping regimes in the 
Bose-Fermi mixture are predicted to be markedly differ- 
ent for repulsive and attractive inter-species interaction 
and highly sensitive to the fermion concentration. For the 
attractive interaction, the growth of the fermionic frac- 
tion has a marginal effect on the bosonic self-trapping, 
but nevertheless leads to the MQST at shorter well sepa- 
rations. This is due to the effective suppression of tunnel- 
ing owing to the narrowing of the bosonic wavefunction. 




FIG. 8: (a) The relative population imbalance between the 
two potential wells in a self trapped state for repulsive 
fermionic component with the concentration r/f, = (solid 
line) and 0.625 (dashed), as a function of well separation, 
(b) Minimum population imbalance in the symmetry-broken 
MQST state as a function of fermionic concentration. Solid 
line is the quadratic fit AN/N To t ~ 1.05 - 3.06r /b - 0.87r 2 fb . 
Inset: the value of well separation corresponding to the min- 
imum imbalance as a function of rft. Solid line is the linear 
fit XMin « 2.88 + I.7I77& 



In the repulsive case, the growth in the fermionic fraction 
has a strong effect on the dynamics of bosons and pushes 
the onset of the MQST regime to larger well separations, 
as a result of the strongly enhanced tunneling due to the 
fermion-induced broadening of the BEC wavefunction. 

Remarkably, the repulsive fermionic impurities can 
cause the population inversion between the ground and 
the first excited states of the BEC in a double- well poten- 
tial and subsequent inversion of symmetry of the MQST 
state. Ultimate consequence of this inversion is the sup- 
pression of self-trapping at large well separations and a 
non-monotonic character of the population imbalance be- 
tween the two wells as a function of well separation. The 
latter is a striking feature of the self-trapping effect in a 
Bose-Fermi mixture that can be detected experimentally. 

Both the symmetry breaking and the suppression of 
self-trapping in the BEC cloud mixed with degenerate 
fermions signal the existence of the new regimes of the 
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dynamics and switching of BECs in atomic waveguides 
and nonlinear interferometers, which can potentially be 
explored in experiments on an atom chip [25j. They are 
also expected to have profound consequence for the for- 
mation and dynamics of the self-trapped states in the 
Bose-Fermi mixtures loaded into periodic potentials. Be- 
yond the mean field, these effects may have an implica- 
tion for the onset of the superfluid to Mott insulator (MI) 
transition in a lattice potential [26|, [2?], [28| , leading, in the 
case of repulsive interaction, to the inhomogeneous sup- 
pression of the MI regime and phase separation. These 
phenomena require further investigation. 
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APPENDIX: VARIATIONAL SOLUTIONS FOR C s 

AND Cox. 

The dynamics of the system Q and the population 
imbalance characterising a self-trapped state (fT0|) are de- 
termined by the interaction dependent coupling param- 
eters, Cj, Coi, and the mode energies, fij. To compute 
them, we use the variational method with the trial func- 
tion given by Eq. (jlip . The mode energies are found from 
the following functional: 



With the trial function (fTTjl . and the fermionic density 
in the form |(SJ) all the components of the Lagrangian 
density can be calculated analytically: 



cir fb aAl 



f da:£/[$o,i] : 
[ dsA[*o,i] = ±%^[(2^-4>o,iT< 1 ] 

JVL za 0,l 



dx £ 6 [$o] 



fj. a Aly/w (1 + w ) 
<7a A^y/ir 

kaoAl^/n 



(l + 3m 2 + 4 W q /2 ) 



2 

,2 , 



4;r 2 erf ( £0 I _ a 2 _ 4a ,2 

a 



(a + 2x ) w H -^w 



f dx 

Jr 



HiaiA\yfK{l - W\) 

aa\A\^JH 



2V2 
ka\A\ypK 



(l + 3wj - Awl /2 



4xperf 



a? 



■4a£ 



+ (a 2 + 2x 2 ) wi] , 



where, ii>o,i = exp(— x\ja^ t ) and the square amplitudes, 
igj = 1/ [2(20,1 i Wo t i)]. The coupling coefficients 
Co, Ci, and Coi can now be computed explicitly as a 
function of the widths ao and a%: 



-2 max 



dx C[<S>j 



H=°. 



(A.l) 



where the Lagrangian density corresponding to the 

equations of motion (O for each nonlinear mode can be 
separated in three parts: C = Ck + Cb + Cf, The ki- 
netic energy contribution {Ck), the bosonic (£&) and the 
fermionic {Ct) mean field contributions are: 



Co = 



Ci 



1 



3m 2 



- 4m 



3/2 



2V27ra (1 + w ) 



2 ■ 



4w 



3/2 



Coi — 



2V27rai (1 - w x ) 

1 - 2w QWl - 2w a *°w 1 - 2WQW1 1 + wl^wj 11 
2^n (a + oT) (1 + w ) (1 - wi) 



C b 



c f = 



2\dx) ' 
l -[^-k{\x\-x f] $ 2 -^> 4 , 



where, 70,1 = a 2 , 1 /(a§ 



The computation of all 



quantities follows from numerically minimizing the func- 
tional (lA.lj) with all numerical constants fixed, and 00,1 
as variational parameters. In the limit of r/f, — > our 
variational results are consistent with the numerical re- 
sults for the pure BEC system presented in [7[. 
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